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b\ 

' We calculate the spin dependent Fermi liquid parameters (FLPs), single particle 

o : 

energies and energy densities of various spin states of polarized quark matter. The 

c : 

expressions for the incompressibility^) and sound velocity (ci) in terms of the spin 
I/"") \ dependent FLPs and polarization parameter (£) are derived. Estimated values of K 

and ci reveal that the equation of state (EOS) of the polarized matter is stiffer than 
the unpolarized one. Finally we investigate the possibility of the spin polarization 



(ferromagnetism) phase transition. 
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00 ■ I. INTRODUCTION 

q ' 

One of the important research areas of the contemporary high energy physics has been the 
study of matter under extreme conditions. Such a matter, in the laboratory can be produced 
by colliding heavy ions at ultra-relativistic energies. Due to asymptotic freedom of quantum 
chromodynamics (QCD), it is predicted that the hadronic matter at high temperature and/or 
density can undergo a series of phase transitions like confinement-deconfinement and/or 
chira. phase transition EL II - In the h,gh density regime QCD predicts the e*,stence of color 
superconducting state |3|, U, |5j. These apart, the possibility of spin polarized quark liquid 
i.e. the existence of ferromagnetic phase in dense quark system has also been suggested 
recently with which we are presently concerned j^, 6]. The properties of dense quark system 
are particularly relevant for the study of various astrophysical phenomenon. 
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The part of the motivation to study the ferromagnetic phase transition in dense quark 
matter (DQM), as mentioned in jf| is provided by the discovery of 'magnetars' |3] where an 
extraordinarily high magnetic field ~ 10 15 G exists it is argued, that the origin 

of such a high magnetic field can be attributed to the existence of spin polarized quark 
matter 9]. To examine the possibility of ferromagnetism in DQM in ref. (| a variational 
calculation is performed where it is observed that there exists a critical density below which 
spin polarized quark matter is energetically favorable than unpolarized state. Subseque ntly 
various other calculations were also performed to investigate this issue 

m 

For example, in |j] it is shown that there is no contradiction between color superconductivity 
and ferromagnetism and both of these phase can co-exist. In lOj, the same problem was 
studied in the large N c and Nf limit while keeping N c /Nf fixed where it was shown that spin 
polarized state can exist, however, in presence of magnetic screening, color superconductivity 
or dense chiral waves disappear. It might be mentioned that such screening is now supported 



by the lattice calculation 



10 



In[5( it is analytically shown that, if quarks are massless, 
ferromagnetism does not appear which is consistent with the conclusion drawn in IOJ. Ref. 8| 
shows that ferromagnetism might appear in quark matter with Goldstone boson current 
where the magnetization is shown to be related to triangle anomalies. 

In the present work, we apply relativistic Fermi liquid theory (RFLT) to study the possi- 
bility of para-ferro phase transition in DQM. The relativistic Fermi liquid theory was devel- 
oped by Baym and Chin where it has been shown how the various physical quantities 
like chemical potential incompressibility (K), sound velocity (ci) etc. can be expressed 
in terms of the Landau parameters (LPs) calculated relativistically. However, the formalism 
developed in [ljj] is valid for unpolarized matter and LPs calculated there are spin averaged. 

In this paper we extend the formalism of RFLT and the required LPs are calculated 
by retaining their explicit spin dependencies. As a result, here various combination of 
parameters like f^i, foi, foi and corresponding; to scattering involving up-up, up- 
down, down- up or down-down spins are appear 14[]. Once determined, these parameters are 
used to calculate quantities like chemical potentials for the spin up and spin down quarks 



or the total energy density of the system as a function of £ 



n )/n q and n q together 



with various other quantities as we shall see. Here n+ and n~ correspond to densities of 
spin up, down quarks respectively and n q = n+ + n~ denotes total quark density {(J. We, 
also compare some of our results with those presented in 6| where more direct approach was 
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adopted to calculate the total energy density from the loop. In addition, the present work 
is extended further to estimate incompressibility and sound velocity in dense quark system 
for a given fraction of spin-up or down quarks. 

Furthermore, in dealing with the massless gluons, we find that naive series expansion fails 
and one has to use hard density loop (HDL) corrected gluon propagator to get the finite 
result for the LPs involving scattering of like spins [ijj]. This however does not cause any 
problem for the calculation of various physical quantities like chemical potential, exchange 
energy, incompressibility etc. We shall see, even though f and f± (suppressing spin indices) 
individually remain divergent, what appears in our case is the particular combination of 
these parameters where such divergences cancel. 

The plan of the paper is as follows. In Sec. II, as mentioned before, we extend the formal- 
ism of RFLT to include explicit spin dependence. In Sec. Ill, we derive spin dependent LPs 
due to one gluon exchange (OGE) for polarized quark matter. Subsequently, we calculate 
chemical potential and energy density. We find the density dependence of incompressibility 
(K) and first sound velocity (ci) with arbitrary spin polarization (£). To compare with 
ref-j^l, we present ultra- relativistic and non-relativistic results and studied para-ferro phase 
transition of quark matter. Sec. IV is devoted to summary and conclusion. In Appendix, we 
calculate various LPs for unlike spin states of scatterer. 



II. FORMALISM 

In FLT total energy density E of an interacting system is the functional of occupation 
number n p of the quasi-particle states of momentum p. The excitation of the system is 
equivalent to the change of occupation number by an amount 5n p . The corresponding 



energy density of the system is given by [13|, LL4j , 

E = E (^S>^+2 2_, J ^y {27r y f P s, P 's'Sn ps Sn P 's', (1) 

where E° is the ground state energy density and s is the spin index, and the quasi-particle 
energy can be written as, 



° P s + ^2j jr^yJps, P 's>5n p , s ,, (2) 
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where e ps is the non-interacting single particle energy. The interaction between quasi- 

particles is given by / pSjP ' s ', which is defined to be the second derivative of the energy 

functional with respect to occupation functions, 

5 2 E 



fi 



ps,p's' 



5n ps 5n p > s > 



(3) 



Since, the quasiparticles are well defined only near the Fermi surface, one assumes 

e ps = n s + v s f (p-p s f ). (4) 

In FLT^theinteraction pa—, W ,„ is expanded on the basis of Legend, po.yno- 
mials, Pi |13l . |14| . The coefficients of this expansion are known as FLPs, which are given 
by 



fr =(2j+i) 



4tt 



Pi{cos9)f } 



pS,p'S' 5 



(5) 



where 9 is the angle between p and p', both taken to be on the Fermi surface, and the 



integration is over all directions of p 



14(. Note that unlike [13j,ll4|, here we retain explicit 



spin indices without performing spin summation. We restrict ourselves for I < 1 i.e. /g and 
f[, since higher I contribution decreases rapidly as the scattering is dominated by the small 
angles and the series converges, here, ff = ~ J2 S > ff s 151 ] . 

The Landau Fermi 
ing amplitude via 
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iquid interaction f pS)P > s ' is related to the two particle forward scatter- 



F F° 
Cp Up/ 



ps,p's' i 



(6) 



where m q is the mass of the quark and the Lorentz invariant matrix Ai pStP > s > consists of 
the usual direct and exchange amplitude, which may, therefore be evaluated by conventional 
Feynman rules. The dimensionless LPs are defined as Ff = N s (0)fi [l^], where iV s (0) is 
the density of states at the Fermi surface is given by, 



N s (0) 



6{e ps - ijl J 



(27T) 



QdegP} ( dp 



2tt 2 



de 



ps / p=pj 



QdegPfSf 
27T 2 ' 



(7) 
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Here gd eg is the degeneracy factor. In our case gd eg = N c Nf where N c and Nf are the color 
and flavor index for quark matter. For spin up (+) and spin down (— ) quark, density of 
states will be change accordingly. In the above expression (dp/de ps ) p=p s is the inverse Fermi 
velocity (1/v related to the FL parameter F(, 



(dp/de ps ) p=p s = (/*7pJ)(1 + F//3). (8) 



v f 



With Eq.([7j) and Eq.® one reads the general relation as 



(9) 



The compression modulus or incompressibility (K) of the system is defined by the second 
derivative of total energy density E with respect to the number density n q , is given by 

0, 12, HQ 20] 



K = 9n„— . (10) 

Now we introduce a polarization parameter _£ by the equations, nl = n q (l + £)/2 and 
n~ = n q (l— 0/2 under the condition < £ < 1 [6]. The Fermi momenta in the spin-polarized 
quark matter then are p^ = + O 1 ^ 3 an d Pj = P/(l — O 1 ^ 3 ) where pj = (vr 2 ?^) 1 / 3 , is 
the Fermi momentum of the unpolarized matter (£ = 0). So, there are two Fermi surfaces 
corresponding to spin-up (+) and spin-down (— ) states, such that E = E(n q ,n~). We have 



dE 

dn a 



dE dn+ dE dn~ 



dn+ dn q 



dn dn q 



- [(l + £)/i + + (l - fa-] 



Using Eq. fflTj) . the incompressibility becomes 



20| 



(11) 



K 



9nq 

4 
9n q 

4 



:i+o 2 l^ + (i-o 2a 



'i + 2 f^ 
1 ^ \n+(o) 



dn 



<i J 



s 



iv-(o) 



(12) 
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where 



l + F* 



dn s q N s (0) 



(13) 



Similarly, the relativistic first sound velocity is given by the first derivative of pressure P 
with respect to energy density E. Since P = 'Yl ls ^ s n s — E 0,0, we have, 



dP 
dE 



dP dn q 
dn q dE 



i + £Kf$ + (i-£K^ 



1 dn~ 



;i + 0^ + + (i-0/^- 



n n 



2[(i+ev + +(i-ev- 



o 2 



N+(0) 



jv-(o) 



.(14) 



In the above Eq. fll2|) and Eq. (|14p . N ± (0) and F^ 1 corresponds to density of states at Fermi 
surface and dimensionless LP for spin up (+) and spin down (— ) quark respectively. For 
unpolarized matter, £ = implying /i + = fi~, Fq = F and JV+(0) = N~(0). From Eq.([T2]) 



and fTTH we have the well known result as K = 9n 



Qdn a 



161 ] and cf 



n q dfi 
/j, dn q 



III. LANDAU PARAMETERS FOR POLARIZED QUARK MATTER 



In this section we calculate LPs for quark matter with explicit spin dependencies. We 
choose spin s along z axis i.e. s = (0,0, ±1) and represent spin-up and down states by 
their signs. For a four- dimensional description of the polarization state, it is convenient to 
define a 4- vector a M which, in the rest frame of each quark, is same as the three-dimensional 
vector s; since s is an axial vector, a M is a 4-pseudovector. This 4-vector is orthogonal to 

(m g ,0)); in any frame we 



the 4-momentum in the rest frame (in which = (0,s),P tl 
therefore have a^P^ = 0,0 ,|22|. 



The components of the 4-vector a M in a frame in which the particle is moving with 
momentum p are found by a Lorentz transformation from the rest frame 22j , 
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with e p = y/p 2 + trig . We can define projection operator P(a) on each of spin polarization, 
7 5 j?£). Accordingly the polarization density matrix p is given by the expression 



P(a) 



1(1 



p(P, s) 



2m, 



iP + m q )P(a), 



(16) 



which is normalized b y t he condition, Trp(P, s) = 1. The mean value of the spin is then 



given by the quantity 
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^^Tr(p 7o S) = ^^Tr(p 757 ) 



2 e n 



p(s-p) 



(17) 



m q (e p + m q ) 

which is reduced to s a 

We consider the color-symmetric forward scattering amplitude of the two quarks around 
the Fermi surface by the OGE interaction. The direct term does not contribute as it involves 



|s in the non-relativistic limit. 



trace of single color matrices like 
from the exchange (Fock) term 6j: 



rA a , which vanishes. Thus the leading contribution comes 



i j 

= l {P ^ pi)2 niAP^P(p',s% (is) 

where a, (3 is the flavor level, i,j is the quark color index, t a (= A tt /2) is the color 
matrix and g is the coupling constant. Since gluon is flavor blind, the u— channel diagrams 
contribute only when a = f3; i.e. scattering of quarks with same flavor 231] . This means 
that the Fermi sphere of each flavor makes an independent contribution. Thus the potential 
energy receives a factor Nf. On the other hand, the quarks with different colors can take part 
in the exchange process, giving rise to a factor N%. Eventually the potential energy density 
is proportional to NfN^g 2 . For the kinetic energy density, there arises an overall factor 
N c Nf. Thus, the factor N c Nf factorizes out of the total energy density and the competition 
between the kinetic and potential energies is not influenced by the number of flavor. The 
number of flavor neither encourages or discourages ferromagnetism 101 ] . 

Without loss of generality, for the calculation of energy density and other related quan- 
tities, we consider one-flavor quark matter. With the help of polarization density matrices 
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given in Eq. ffTB]) . we have from Eq. fllSI) the interaction amplitude as 
2n 2 1 

x{m 9 (£:p + m,)(p' ■ s)(p' • s') + m q (s p , + m q ){p • s)(p • s') + (p ■ p')(p ■ s)(p' ■ s')]. 

(19) 

From Eq.Q the quasiparticle interaction parameter is given by 

/■ex '^q \ azx 

ps,p's' = ~~ JV[ ps,p's' 
£p £p' 

Here the spin may be either parallel (s = s') or anti-parallel (s = — s^. Thus scattering 



possibilities are denoted by (+,+), (+,—), (— , — ) etc. Motivated by [15|, in analogy with 



isospin we define spin dependent interaction parameter as f pp , = \{f pp t + fp p i ) and f ppl = 

UfpV + Note that > fp P > = fp7- 

For (+, +) scattering the interaction parameter is given by 

f + t , + — To~ TsT, 77 [2m„ — p| 2 (l — COS#) — pt 2 COS #i COS #2 

j pp'\p=p'=p) 9 £ + 2 p+ 2 (l -cos#) L q f V ' f 

+ , t — {mJet + rn q )pt 2 (cos 2 6>i + cos 2 2 ) +p| 4 cos0cos#i cos# 2 }], 

(ej+m g ) 2 1 1 1 

(21) 

where p ■ s = cos#i ; p' • s — cos6 2 and Fermi energy = (p^ 2 + m 2 ) 1//2 . Since spin and 
momentum have no preferred direction, we have done angular average of the spin dependent 



parameter 241 ] : 



pp i\p=p'=p+ 



47T 



; /•++ 



2 



47T J pp'\p=p'=pj 
2 

'2m 2 



9£j 2 pJ 2 (l - cos I 



COS I 



2m q p J f 



3(5" 



(22) 



1 With the help of Eq. (jSJ) along with the Eq. (1221) one can find LPs, but it is to be noted that 
fo~i or f~~ are individually divergent because of the term, (1 — cos 8), in the denominator of 
the interaction parameter. This divergence disappear if one uses Debye screening mass for 



[1] denoted hereafter f pp > = f 1 



pp' ■ 
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gluons or e quiv alent ly use HDL corrected gluon propagator while evaluating the scattering 
24j. Note that the combination ( / ++( ) 



amplitudes 



12 



1 f ++(— ) 



is, however, finite as 



in this case the divergences cancel and we do not calculate the LPs separately. It would, 
however, be interesting to see how do the results modify if HDL calculations are performed to 
evaluate f^^ , f + [' and the corresponding physical quantities. The numerical estimates 
suggest that for the results what we present here, the effect of HDL corrections are expected 
to be small. 
From Eq.([SD, 



f++ _ 1 f + 

JO gil 



+ 1 
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f Pf J -i 

2 



2mt 



COS I 



2m q p^ 2 
3(e^ + m q 



d(cos i 



2m\ - p + f 2 + 



2m q pY 
3(£j + m q 



(23) 



The above combination will appear in the calculation of the chemical potential and other 
relevant quantities. For (+, — ) scattering, the angular averaged interaction parameter yields 



p=pj,p=p f 



9epj 



+2 

m q pj 



rn 



+ 



3(V^ + m q ) 3(e f + 




(m q — EfSf + PfPf cos ( 



(24) 



It is to be noted that, individual LPs for scattering of unlike spin states are finite i.e. free 
of divergences, in contrast to the case involving scattering of like spin states (For details see 
Appendix). 

A. Chemical potential 



Now we proceed to calculate chemical potential, which, in principle, will be different for 
spin- up and spin-down quarks, denoted by fi s with s (or s') = +, — for matter containing 
unequal densities of up and down quarks. To determine the chemical potential with arbitrary 
polarization £, we take the distribution function with explicit spin index (s or s'), so that 



variation of distribution function gives 



13 



20 



5n" q 



-AT s (0) 



(25) 
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where N s (0) is given by the Eq.([7j). The Eq. (l2Sl) yields 



dni N s (0) V dn 

H S 

Separately for spin-up and spin-down states we have 



" 1 -E^- (») 





iV-(O) 



/o" 



(27) 



where the superscripts ++ and H — denote scattering of quasiparticle with up-up and up- 
down spin states. For unpolarized matter the upper and lower component become equal 
which gives rise to the well known result fl^ | 



/id/i 



. 9degfJ>Pf, 



d Pf . 



(28) 



In general the chemical potential (both for spin-up and spin-down) is the combination of like 



and unlike spin states. By adjusting the constant of integration 
of spin-up quark turns out to be 



14J, the chemical potential 



6lT 2 E~, 



11 2, 

— m q \n 



Pf 



m n 



pjm q - 



72tt%" 



2m q 
p+ 



In 



pj+Pf 
P + f~P~f 



-14m?ln ^ 



2m, 



Pf 



Pf 



3m n pf In 



In 



Pf £ f 
2 

P~f + Pf 



In 



-%2m; + 3 P ; a )ln T 7+ g^ ) 
Pf \Pf £ .f Pf e f 



Pf ~Pf 

Pf +Pf 
p+-pj 

Pf 



Pf £ f +Pf £ f 
Pf e f-Pf £ f 



Pf 



p f 2 }\n 



£ f £ f 



m' 



~PfPf~ \ 



£ f £ f 



mi 



PfPf) 



6p f e 



(29) 



In the above equation the term in the first square bracket arises due to the scattering of 
like spin states (++), while the latter comes from the scattering of unlike spin states (H — ). 

Similarly, for spin-down quark, one may determine \x~ by replacing p^ with pj and Ej 
with sj in Eq.(|25jl. 



± 



11 

For the numerical estimation of the above mentioned quantities, following ref . {(], , we 
take a c = g 2 /4:7i = 2.2, is the fine structure constant of QCD and m q = 300MeV. In 
Fig([T]) we plot chemical potential for spin-up and spin-down quark as a function of density 
with order parameter £ = 0.5. In real astrophysical calculations, the chemical potentials 
are determined by the /3-equilibrium conditions where the condition of charge neutrality is 
also imposed. In FigflTJ), we, however, use density n q and polarization parameter £ as input 
parameters and Eq. (l2U|) is used to determine \i for a system with one flavor. 
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FIG. 1: Density dependence of chemical potential of spin-up and spin-down quark denoted by solid 
and dashed curve respectively. 



B. Energy density 



Once the ji is determined, one can readily calculate the exchange energy density by 



evaluating [14, 
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E ex = J dn q ((i-£f) 



(30) 



After summing up over the color degrees of freedom and evaluating over the Fermi sur- 
faces, we have the exchange energy density. The latter consisting of all types of scattering 
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amplitudes, can be written as 



E P 



E++ + E-- + e: 



(31) 



which we evaluate numerically. The total kinetic energy density for spin-up and spin-down 
quark is given by 



E 



km 



167T 



s=± 



Pf £ f( £ f+p s f)- m i ln 



(32) 



where e s j- = (p s j + m 2 ,) 1 / 2 . The total energy is given by the sum of the kinetic energy and 
the interaction energy E ex i.e. 



E 



tot 



Ekin E ex . 



(33) 



Now we calculate incompressibility and sound velocity by using Eq. f)12l) and Eq. fll4l) . In 
Fig.([2]) and Fig.Q we plot the density dependence of incompressibility and sound velocity. 
This shows for higher value of the order parameter £, the incompressibility and the sound 
velocity becomes higher for the same value of density. Thus the EOS for polarized quark 
matter is found to be stiffer than the unpolarized one. 




FIG. 2: Incompressibility K in quark matter as a function of density for different polarization 
parameter. 
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FIG. 3: First sound velocity ci in quark matter as a function of density for different polarization 
parameter. 



C. Phase transition 



Bloch first pointed out the possibility of ferromagnetism o 



exchange interaction induces spontaneous spin polarization 27| . Consider the spin polarized 



electron gas where the Fock 



electron gas interacting by the Coulomb interaction in the background of the positively 
charged ions. Since the direct interaction gives no contribution due to charge neutrality, 
the Fock exchange interaction gives the leading contribution as the interaction energy. For 
spontaneous ferromagnetism, the interaction energy dominates over the kinetic energy 6|, 



21 



28]. 



Therefore, if the exchange energy due to OGE interaction is negative and becomes greater 
than the kinetic energy at some density, the quark matter becomes polarized giving rise to 



ferromagnetism 



6]. 



To check whether our results for the total energy density are consistent with ref. we 
consider two limiting cases corresponding to the ultra-relativistic (UR) and non-relativistic 
(NR) regimes. In the ultra-relativistic (UR) limit, pj ^> m q , then using Eq. fl29l) we have 

/i+' ur = p) 



3vr 



Pf 



P? 



Pf 



(34) 



Similarly one can find /i ' ur by replacing p^ with pj. 
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One can arrive at the same expression fi ±,ur by taking UR limit of the scattering ampli- 
tude. For (+, +) scattering one gets the interaction parameter as 



H — h,ur 



J pp' 



9pp' 



'1 + cos 6*i cos 



(35) 



After taking angular average of the interaction parameter and with the help of Eq.((5]) 
we find that f± + vanishes. Thus we have 



++,ur 



J pp 1 



H — \-,ur 




9pf 



Similarly for 



p=p'=pj ciyj 

scattering, the interaction parameter yields 

-.2 

'1 — COS#i COS 9 2 ) ■ 



H — ,ur 



J pp' 



g- 



9pp' 



(36) 



(37) 



The only existing LP is f$ and other higher order LPs does not contribute. Hence we get 

,2 



r + - 
J pp' 



9 



(38) 



p= P j,P'=p f ^PfPj 
It is observed that, in UR limit, all the LPs are finite. Now the chemical potential for spin-up 
quark is found to be 

-2" 



/i+' ur 



a,. 



Pf + — 



3tt 



Pf 



P_l_ 

P + f 



(39) 



The chemical potential, \x ' UT \ can be obtained by replacing p^ with pj in Eq. (l39l . 
Using Eqs. (j3Dj) and (J3*T]) . the exchange energy densities are given by 



Tp++,ur 
ex 



i^Pf 



771 ,ur _ Clc j,- 1 

^ex — 8w*Pf 



E+ 



ur _ _ac_ „+2 -2 
- 4w 3Pf Pf ■ 



(40) 



Thus the final expression for the exchange energy density in the UR limit is found to be 



El 



a, 



8tt 



■pj [(i+o 4/3 +(i-o 4/3 +2(i-a 2/3 ]- 



(41) 



This result is same as in ref.jfi]. 

Similarly from Eq. (1321) . the kinetic energy density in UR limit takes the following form 

3: 



rpur 



3p 



f 



8vr 2 



[(i + 4/3 + (i-0 4/3 ] 



(42) 
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In the NR limit, p(or p') <C m q , the interaction parameter reduces to a simple form, 



J ps,p's' 



9pp' 



1 + 8 -a? 

[1 - cos#) 



(43) 



For spin anti-parallel interaction s = —s', then fpg p i s i = 0. Thus the contribution due to 
the scattering of quarks with unlike spin states vanishes and the dominant contribution to 
energy density comes from the parallel spin states (s = s'). For the (s, s) scattering, the 
interaction parameter yields 

■p'=P s f 9pf(l - cos 9) ' 



rnr, 
J pp' 



(44) 



where s = + or — according to scattering process. In NR limit one gets 



(/o 3^ 1 ^ 



The NR chemical potential fi nr is given by 



2g 2 



9p s f 



(45) 



m n 



g- 



q 3tt 2J 

Using Eq. fl30|l . the exchange energy density for the (+, +) scattering is given by 



(46) 



Similarly for 



scattering, we have 



El 



9 4 



8tt' 



8tt 



p 4 /(i + 4/3 - 



//(l-0 4/3 - 



(47) 



(48) 



As in the NR limit E+ x = E ex *~ = as mentioned before, so from Eq. (l3~T|) the exchange 
energy density yields 



jpnr 



a. 



2tt 3 



rf[(i + 4/3 + (i-0 4/3 ] 



(49) 



Thus the energy density, in this limit, becomes negative. The kinetic energy density turns 
out to be jfl 

,5 



rpnr 



207r 2 m o 



[(i + 5/3 + (i-0 5/3 ]- 



(50) 



In NR limit ferromagnetism can appear as a consequence of competition between the kinetic 
energy and the Coulomb potential energy 10|. The latter favors the spin alignment due 
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to quantum effect. When the energy gain due to the spin alignment dominate over the 
increase in the kinetic energy at some density, the unpolarized state suddenly turns into 
the completely polarized state 29(. On the other hand, in UR limit, the contribution to 
the energy density not only comes from the like spin states but also unlike spin states of 
scatterer (see Q| for detailed discussion). 
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FIG. 4: The total energy of quark liquid as a function of polarization parameter at n q = 0.1 fm 



0.14 fm 3 and n„ = 0.2 fm 3 . The critical density is found to be n c g = 0.14 fm 3 in this case. 



To check the consistency we compare our result derived in RFLT approach with that 
of ref 6J derived from two loop ring diagrams. In Fig.(|3J) we plot E tot /n q as a function of 
polarization parameter £. The results clearly show that for lower density (< 0.14 fm~ 3 ), 
total energy favors at £ = 1 which indicates completely polarized state; while at higher 
density, the system becomes unpolarized (£ = 0). Thus the polarization parameter suddenly 
changes from £ = 1 to £ = as one increases the number density of the system. So the 



phase transition is first order and the critical density n c q is around 0.14 fm 



-3 



In Fig. ([5]) we show total energy as a function of polarization parameter for different 
densities. In every plot, there is a minima which corresponds to a possible metastable state. 
We notice that when density increases metastable state arises for lower values of polarization 
parameter £. For example, at density ~ 0.2 fm -3 minima arises at £ = 0.1 while at density 
~ 0.35 fm -3 minima arises at £ = 0.03. Thus the metastable state shows a tendency of 
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(a) (b) 




(c) (d) 



FIG. 5: Metastable ferromagnetic state as a function of polarization parameter for different densi- 
ties. 

disappearance as the density increased. 

IV. SUMMARY AND CONCLUSION 

To summarize and conclude, in this work we have applied RFLT to study the properties 
of dense quark matter. Accordingly, we calculate the FLPs by retaining their explicit spin 
dependencies. We also show how the physical quantities like chemical potential of spin 
up and spin down states, their energy densities and the quantities like incompressibility, 
sound velocity for polarized quark matter can be expressed in terms of these spin dependent 
RFLPs. For the scattering involving like spin states, the LPs f$i and are found to 
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diverge. However, we show that the combination in which they appear in the calculation of 
the physical quantities such divergences cancel. For the scattering involving unlike spin states 
no such divergence appear. The appearance of such divergences is related to the unscreened 
gluonic interaction between the quarks which might be cured by invoking hard dense loop 
corrected gluon propagator. We do not perform such calculation here and postpone this for 
future investigation. As far as the equation of state (EOS) is concerned, we in the present 
model find that the EOS for the polarized quark matter is stiffer than the unpolarized one. 
In addition, we also show that there exists a metastable state which disappear at higher 
density, although it seems that the effect is tiny. 

We reconfirm that DQM can exhibit ferromagnetism at low density as was originally 
suggested in 6|. However, the density at which the spin polarized ferromagnetic state in 
the present model might appear depends strongly on the quark mass. The critical density 
increases with increasing mass. In Fig. (J3J), we observe that states with £ appear only 
below or around normal nuclear density where deconfined quark matter is not likely to 
exist. We cannot, however, ascertain the critical density from the present analysis where we 
restrict ourselves only to OGE diagrams and one flavor system. In this regime, multi-gluon 



exchange processes [10| might play an important role. Furthermore, the correlations as given 
by the ring diagrams can also change the conclusion. Further work therefore is necessary to 
understand the existence of ferromagnetic quark matter in real multi flavor system which 
might appear in astrophysics. 

V. APPENDIX 

In the text the interaction parameter f^ p 7 for unlike spin states was calculated. Here we 
give detail expression of Landau parameters. With the help of Eq.((5]), the LPs are given by, 



Jo 



lSe+eJ 



3pjp f (m q + ej)(m q + e f ) ' \m 2 q +p]p f -e]s f 



(51) 



and 
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ft 



m q [m q 



_ 2m g [e / p+ + m q (pj + p f )+ejp f 
p^pj(m q + e^)(m q + ej) 

- e+ej) [ejpf + m q {pf + pf) + 



P/ P f (m q + ej)(m q + e f ) 



Using Eq. (l5"Tj) and Eq. (152ft we have 



In 



m g +PfPf 



— eje 



PfPf 



£ f £ f 



(52) 



/o g fi 



+2 

m q pj 
3(e^ + m q 



PfPf 



+ 



711 



+ 



3(e f + m q 



x 



(p7V) 



in 



P/P/ 



• (53) 
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